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The s-wave interaction is usually the dominant form of interactions in atomic Bose-Einstein con- 
densates (BECs). Recently, Feshbach resonances have been employed to reduce the strength of 
the s-wave interaction in many atomic speicies. This opens the possibilities to study magnetic 
dipole-dipole interactions (MDDI) in BECs, where the novel physics resulting from long-range and 
anisotropic dipolar interactions can be explored. Using a variational method, we study the effect of 
MDDI on the statics and dynamics of atomic BECs with tunable s-wave interactions. We bench- 
mark our calculation against previously observed MDDI effects in 52 Cr with excellent agreement, 
and predict new effects that should be promising to observe experimentally. A parameter of mag- 
netic Feshbach resonances, tdd.max, is used to quantitatively indicate the feasibility of experimentally 
observing MDDI effects in different atomic species. We find that strong MDDI effects should be 
observable in both in-trap and time-of-flight behaviors for the alkali BECs of 7 Li, 39 K, and 133 Cs. 
Our results provide a helpful guide for experimentalists to realize and study atomic dipolar quantum 
gases. 

PACS numbers: 03.75.-b, 67.85.Bc, 67.85.De 



I. INTRODUCTION 

The physics of ultracold dipolar quantum gases is a rich 
and promising area of research. There is great interest 
in the physical behaviors that result from dipole-dipole 
interactions [TJ [2]. An ultracold atomic gas of atoms 
that posses a magnetic moment will have a magnetic 
dipole-dipole interaction (MDDI). However, it is often 
difficult to observe the effects of MDDI in many atomic 
species (particularly alkalis), where the MDDI is typi- 
cally much weaker than the isotropic, s-wave interactions. 
For atomic species with magnetic Feshbach resonances, 
however, the s-wave scattering length, a s , can be tuned 
via magnetic fields [3J. Employing Feshbach resonances 
has led to fruitful and impressive developments in ultra- 
cold atom research. Of interest here is that it allows for 
the exploration of MDDI in Bose-Einstein condensates 
(BECs). By tuning a s to near zero, MDDI can become 
the strongest interaction in the BEC. 

The effect of MDDI in ultracold atomic gases was first 
observed using a BEC of 52 Cr atoms, which posses a 
strong magnetic moment, fi, of 6 Bohr magnetons (/xs) 
[3]. The MDDI were observed to affect the aspect ratio 
of the 52 Cr BEC in time-of-flight (TOF) free expansion. 
The MDDI effect on the stability of a BEC was also 
experimentally studied in 52 Cr, where for various trap 
configurations the MDDI either made the BEC more or 
less stable against collapse [5]. More recently, in 7 Li the 
MDDI effect was seen when comparing the axial length 
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of the BEC near a s — for two different trapping ge- 
ometries [6]. Effects of MDDI have also been observed 
on the decoherence rate in a 39 K BEC atomic interferom- 
eter [7], and on the spin domains of a spinor 87 Rb Bose 
gas [5]. Very recently, a BEC has been realized with 
164 Dy, with a strong magnetic moment of IQ^b, exhibit- 
ing dipolar effects even with no tuning of the scattering 
length [9 . Many further effects due to MDDI have been 
predicted such as the excitation of collective modes by 
tuning the dipolar interaction |10j . the emergence of a 
biconcave structure with local collapse [IT] [T3J, and the 
modification of the phase diagram of dipolar spin-1 BECs 
[B] [13] [IS], of the soliton stability in ID BECs [H], and 
of vortices in BECs [17] . We refer the reader to recent re- 
views for a further discussion of the multiplicity of MDDI 
effects [JJ [SJ . 

Motivated by such rich physics, we theoretically model 
the effects of MDDI and possibility of experimentally de- 
tecting such effects in BECs of 52 Cr and all the alkalis. In 
Section II, we explain the variational method we employ 
to model MDDI using a cylindrically-symmetric, Gaus- 
sian Ansatz for the BEC wave-function [T5]. In Section 
III, we start by discussing the relevant parameters for 
our simulations. We also introduce a key quantity used 
in this paper, the ratio of s-wave scattering length, a s , to 
a length defined for MDDI, add- Next in that section we 
present simulations of the effects of MDDI in 52 Cr, in- 
cluding those reported in Refs. [5] and [19], and predict 
several additional effects that should be readily observ- 
able. We also present simulations for the effects of MDDI 
in 7 Li, 39 K and 133 Cs, the three alkali BECs we identified 
as promising for observing MDDI effects. In Section IV, 
we conclude. 
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II. VARIATIONAL METHOD 

While a few other methods exist to model dipole-dipole 
interaction effects in BECs 1, 2(H23|. the variational 
method we employ has shown great utility because its 
simple, analytic solutions are valid over a wide range of 
experimental parameters [HJ [24] ES] . Two types of inter- 
actions are considered. The s-wave interaction is charac- 
terized by the s-wave scattering length, a S} and for dipole- 
dipole interactions a parameter defined as add is used. 
MDDI can be characterized by add = Mo A* 

M/(12tt/i 2 ), 

where fi is the magnetic moment of the atom, M its mass, 
and /iq is the permeability of free space [55] . The atom- 
atom interaction potential thus has two terms, one for 
each interaction: 
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Using that interaction term, the Gross-Pitaevskii equa- 
tion for a BEC takes the form (in dimensionless units): 
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where the length unit is a r = ^h/mu) r , the trap frequen- 
cies are cu r and w 2 for the respective radial and axial di- 
rection, the time unit t — 2n/u> r , ip is normalized to unity 
(|Vf = 1), and V ext {r) = \ (x 2 + y 2 + Xz 2 )/2 is the trap 
potential where the trap aspect ratio is A = uj z /uj r . For 
this paper, we adopt the nomenclature for trap shapes 
that is common to experiments: the symmetrical a.k.a. 
spherical (A = 1), the cigar a.k.a. prolate (A < 1), and 
the pancake a.k.a. oblate (A > 1) shapes [5]; we also 
assume that the magnetic field is applied along the ax- 
ial direction. Using a cylindrical-symmetric, Gaussian 
Ansatz for the BEC wave-function 1271 
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the variational method results in two differential equa- 
tions that describe the mean axial, q x , and radial, q ri 
lengths of the BEC (detailed solution in Ref. [IB] , noting 
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These coupled differential equations model BEC be- 
havior with (keeping add) and without (setting add = 0) 
MDDI effects, and they can be numerically solved to 
model three experimentally relevant situations: 

1. The static, in-situ sizes for a trapped BEC are 
found by setting the time-dependent components, 
q r and q z , to zero. 

2. In-trap dynamics are modeled by keeping all terms. 

3. Time-of-flight (TOF) free expansion behavior is 
modeled by removing the the terms q r and X 2 q z , 
which represent the trapping potential, on the left 
sides of Eqns. [4] 



III. RESULTS 

In this section, we employ the above method to solve 
for both the in-trap and TOF behaviors of BECs with 
MDDI [29]. We also find the threshold a 3 , denoted 
a threshoi<\ below which the BEC is unstable and collapses. 
We benchmark our variational calculations against avail- 
able experimental results in 52 Cr and find good agree- 
ment (Section III.B). We also make predictions of MDDI 
effects in the alkalis, and find that 7 Li, 39 K, and 133 Cs are 
the species most favorable for the exploration of MDDI 
effects. 



A. Parameters 

The input parameters used in our simulation are the 
number of atoms in the trap (iV), the magnetic dipole 
moment of the atom (/j,), the mass of the atom (M), 
the axial (f z ) and radial (f r ) frequencies of the trap 
(27r/ ri2 = w r ,z)j and the s-wave scattering length (a s ). 
An applied magnetic field can tune a s via a Feshbach res- 
onance with an analytic approximation given by a s (B) = 

o-bg fl — B A B J , where A is the width and is the 

location of the Feshbach resonance, and ab g is the scat- 
tering length far from any resonances. An experimen- 
tal limit for reaching small a s is the precision of control 
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Species 


\F,m F > fi/nJ^ 


a a ,min («o) 


add (a ) 


£dd,max 


'Li 


|1, +1 > 


0.84 


0.0007 


0.041 


58.4 


23 Na 


|1,+1 > 


0.50 


0.572 


0.047 


0.08 


39 K 


|1,+1 > 


0.95 


0.0022 


0.287 


130 


4i K 


|1,-1 > 


0.08 


0.104 


0.002 


.020 


85 Rb 


|2, -2 > 


-0.57 


0.067 


0.223 


3.33 


s7 Rb 


|1,+1 > 


0.73 


5.98 


0.374 


0.062 


133 Cs 


|3, -3 > 


-0.75 


0.0096 


0.607 


63.2 


52 Cr 


3, -3 > 


6 


0.473 


15.2 


32.2 



a Calculated at Bq 

TABLE I. Parameters for variational computations, with the 
maximum value of e^d as calculated from Eqn. [6] 



over the magnetic fields. In typical ultracold atom ex- 
periments an experimental precision SB/B of approxi- 
mately 10~ 5 can be realized [50H52"] . As n depends on 
the strength of the magnetic field, we calculate /i at 
Bq = Boo + A (where a s = 0), denoted Across- For ref- 
erence, some parameters for known Feshbach resonances 
from the literature are listed along with the calculated 
values of Across in Appendix [X] Table |III| 

To compare the potential for experimentally observing 
MDDI effects in the various atomic species of interest, we 
employ the dimensionless ratio €dd = &dd/ a s used in Ref. 
[T5] . focusing on its maximal value that can be achieved 
experimentally: 
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where a Siin i n ~ a b g (SB/ A) is the minimal a s that can be 
achieved given a typical experimental magnetic field sta- 
bility (assumed to be 5B/B w 10~ 5 ). The alkali species 
that are the best candidates for observing and studying 
MDDI effects in BECs are clearly 7 Li and 39 K and 133 Cs 
(see Table [i]). We note that 52 Cr, while having a much 
larger a^d, does not have a broad Feshbach resonance 
that allows for as precise tuning of a s as in some of the 
alkalis. 

As an initial verification of our model, we compute k, 
in-situ (t — 0) and in the TOF asymptotic limit (t ~ > 
oo ) assuming only s-wave interactions, for both strongly 
interaction (a s — > Large) and non-interacting (a s = 0) 
cases. The expected in-trap and TOF behaviors of a 
BEC in such limiting cases are (see Refs. [551 134] ) : 



1. If a s = and t — 0, expect k 

2. If a s = Large, t = 0, expect k = ^f- = A. 

3. If a s = 0, t -> oo, expect k = = A" 1 / 2 . 

4. If A <C 1 or A 3> 1, and if a s = Large with t — > oo, 



expect k ■■ 



2\- 1 /tt. 



Our variational calculation (performed with ji = in 
these cases) does reproduce these expected values for k. 
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FIG. 1. The calculated in-trap BEC aspect ratio (a), radial 
length (b), and axial length (c) of a 52 Cr BEC with (solid, 
blue) and without (dashed, black) MDDI for a range of trap 
aspect ratios, A. In this simulation, the f avg = 700 Hz, the 
atom number N — 2 x 10 4 , and a s = 15ao. The magnetic 
field is, as for all simulations in this paper, aligned along the 
axial direction. 



B. Effects of MDDI in 52 Cr BEC 

The first observation of MDDI in a BEC was with 52 Cr. 
As 52 Cr is the most studied atomic species so far for 
MDDI effects in BECs, we present and benchmark our 
calculation and results for 52 Cr before discussing the al- 
kalis. We used 52 Cr to demonstrate four characteristic 
MDDI effects, discussed in detail below. 



1 . Effect of MDDI on m-situ aspect ratio of a BEC 

A calculated result of a 52 Cr BEC trapped in a har- 
monic trap is shown in Fig. [T] For a nearly symmetrical 
trap (A ~ 1, which is the case for the experiment in 
Ref. [S]) the MDDI increase the axial length of the BEC 
and reduce its radial length compared to a BEC with no 
MDDI, leading to a decreased aspect ratio. However, if 
the trap is very prolate, A <C 1 (very oblate, A ^> 1), we 
find that the BEC will shrink (expand) in both the axial 
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and radial directions, in a way that leads to an increased 
aspect ratio. There are two values of A, one for when 
the trap is oblate and one for when the trap is prolate, 
where the MDDI do not change the aspect ratio of the 
BEC. Similar in-situ effects of MDDI on n were found in 
all simulations of stable BECs for the alkali as well. 
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FIG. 2. Stability diagram and aspect ratio of a Cr BEC, 
with parameters (f a v g = 700 Hz, and N = 2 x 10 4 ) cho- 
sen to resemble those in the experiment Ref. [5] . The BEC 
aspect ratio, k, solved using our method, is plotted in the 
color map as functions of a 3 and A in the stable regime (log 
10 scale). The energy variational solution and experimental 
data from Ref. [5] for a t hreahold are included for comparison. 
This shows the effectiveness of our method in solving not just 
the al hrcBhold where the BEC collapses but also the BEC size 
in the stable regime. 



2. Effect of MDDI on stability of a trapped BEC 

For BECs with only s-wave interactions, it is shown 
that if a s < (attractive interactions) the BEC will col- 
lapse if N > 0.55aho/\a s \ : where ay lo = ^Jh/raui and 
a) is the average trap frequency [3"3H3"rJ] . For any purely 
s-wave interacting BECs held in a three dimensional har- 
monic trap, the BEC will not collapse if a s > 0. However, 
MDDI can destabilize or stabilize a BEC that would oth- 
erwise be stable or unstable. The effect of MDDI on the 
stability of a trapped BEC was first observed in 52 Cr [5]. 
For cigar traps with the B-field aligned along the axial di- 
rection, MDDI can lead to BEC collapse at a larger value 
of a s (in some cases, even a s > 0) than in an otherwise 
identical BEC with no MDDI. For pancake traps, the ef- 
fect is opposite, and the MDDI stabilize the BEC. Koch 
et al. performed an experimental study of this MDDI ef- 
fect [5]. To model their results, Koch et al. employed an 
energy argument based on a variational Gaussian Ansatz 
for the BEC density distribution. Their method relied 
on finding when a minimum in the energy vanishes as 



a s is reduced. Different from their method, our varia- 
tional calculation solves Eqns. [4a] and [4b] over a range 
of a s , and find the threshold a s when those equations do 
not have stable numerical solutions, which indicates BEC 
collapse. Our method allows not only the determination 
of the threshold a s for collapse (which agrees with those 
obtained in Ref. [5]) but also the axial and radial BEC 
lengths over the entire range of a s where the BEC is sta- 
ble. We employed our method, using N and f avg identical 
to those used in Fig. 3 of Ref. [5 , to solve for the thresh- 
old a s and also k for a s > a t hreshold ove r a range of A (see 
Fig. [2| . We find excellent agreement with both Koch et 
aVs calculation and their experimental data. Our calcu- 
lation shows, as observed in the original experiment [5], 
that the anisotropic dipole-dipole interactions cause the 
threshold a s to depend strongly on the trap geometry. A 
similar dependence of a t hreshold n A is also seen later in 
this paper for alkali BECs. 



3. MDDI effect in time- of -flight behavior 

MDDI affects the behavior of a BEC released from a 
trap, and such an effect was first observed in an exper- 
iment on 52 Cr, where the BECs aspect ratio in TOF 
changed when the applied magnetic field was along the 
axial verses radial direction [4] [19]. 

We have simulated the experiment in Ref. [19] (specifi- 
cally the data shown in their Fig. 4) and find good agree- 
ment. We assume cylindrically-symmetric (f r — 600 Hz 
and f z — 370 Hz) trap for ease of calculation with pa- 
rameters approximating their "trap 2" , which had fre- 
quencies ^,^,^=660,540,370 Hz, and N= 3 x 10 4 . The 
values of a s (112 ao, 96 do, 30.5ao, 20.5eto) and resulting 
edd are chosen to match the values in Fig. 4 a-d of [15] 
(e dd = 0.14, 0.16, 0.5, 0.75). Even with the cylindrical 
trap approximation, the simulation results agree with the 
observed data quite well, and clearly shows the effect of 
MDDI to reduce the aspect ratio of the BEC in TOF 
from this trap configuration, as shown in Fig. [3j We fur- 
ther show that in the case of small a s , MDDI- induced 
collapse could also be observed after a BEC is released 
from a trap if the a s is tuned to a small value upon re- 
lease. Fig. [3] also shows the MDDI-induced near-collapse 
(e) and the collapse (f ) of atoms in free-expansion. 



4- Effect of MDDI on m-trap dynamics 

Our model can also solve the in-trap dynamics of BECs 
with MDDI. A simulation of the aspect ratio, radial size 
(g,.), and axial size (q a ) of a trapped 52 Cr BEC — initially 
perturbed from its static state — evolving with time is 
shown in Fig. [4j revealing an oscillatory behavior. One 
effect of the MDDI in this situation is to reduce ampli- 
tude of the oscillations. The oscillations may be difficult 
to observe in-situ due to the small condensate size, but 
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FIG. 3. The aspect ratio of a 52 Cr BEC vs the time-of-flight duration,fTOF, for different scattering lengths (a-d) with parameters 
chosen to resemble those in Ref. [19]. We additionally show simulated TOF behavior where MDDI nearly (e) and actually (f) 
causes BEC collapse. Upon the release of the BEC, the scattering length is tuned to a s — 10 ao in (e), and a s = 5ao in (f). 
For (a-f), N = 3 x 10 4 , f r = 600 Hz, and f x = 370 Hz. 



would become easier to observe in TOF measurements 
taken from different instants of the oscillations. As seen 
in Fig. [4] (d) and (e), the aspect ratio in TOF changes 
by a factor of about 2 because of the MDDI. An in-depth 
treatment of the modes of oscillatory, in-trap dynamics 
of a BEC with MDDI can be found in Refs. [57] and 

MDDI effects have been observed in collective oscilla- 
tions of 52 Cr BECs [31]. While the results of Ref. [35] 
cannot be directly simulated by our method because their 
magnetic field was not aligned on the axial direction of 
the BEC, we are motivated by this experiment to study 
the collective oscillations of BECs with MDDI. By tak- 
ing the Fourier transform of the time-dependent, in-trap 
oscillations, we obtain the frequency spectra of the col- 
lective oscillations of a BEC with MDDI, exemplified in 
Fig [5j The three lowest modes of oscillation were ob- 
served, and shifts due to MDDI were seen in both the 
amplitude and frequency of the modes. 



C. Results for alkalis 

1. Highlight of MDDI effects in Alkali BECs: stability 

A central point of our paper is that MDDI effects are 
also possible to observe in BECs of the alkalis, even with 
a much smaller \jl than 52 Cr. To show this, we compare 



dipolar collapse for BECs of various species in Fig. [6] If 
the MDDI are not included in the model, then the BECs 
are stable for any positive value of a s . However, with 
MDDI the BEC aspect ratio decreases more substantially 
as a s is reduced towards zero. The BEC collapses be- 
yond a a t hreshold ; indicated by the heavy dot. The value 
of a t hrcshold C an be compared with the a s , m i n of Table [ij 
indicated by the colored bars in the figure. This compar- 
ison clearly indicates whether the observation of collapse 
due to MDDI is feasible with current experimental abil- 
ities and reveals 7 Li, 39 K, and 133 Cs as promising alkali 
species for such an observation. 



We find that 7 Li, 39 K, and 133 Cs possess the greatest 
potential among alkali BECs for exploring MDDI effects. 
Several examples of such effects are shown in the sim- 
ulations below. We describe 7 Li in detail, and similar 
results are presented for 39 K and 133 Cs. The |1, +1 > 
state of 7 Li is excellent for studying MDDI effects as it 
has the the widest known Feshbach resonance of the alkali 
atoms; the slope at zero-crossing is only « 0.1an/G [BJ. 
Using our variational method, we find that each of the 
MDDI effects discussed so far in 52 Cr should be observ- 
able within current experimental capabilities with 7 Li as 
well. 
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7 Li BECs [40], and here could be used to reveal MDDI ef- 
fects, (d) shows the in-trap evolution of a BEC initially 
at a s — O.Olao and then tuned to O.OOlao at t = to 
induce in-trap collapse due to MDDI. Neglecting MDDI 
would result in a stable, oscillating BEC, as seen in Fig. [7] 
(d). The results in Fig. goffer four methods for detecting 
MDDI in 7 Li BECs. Similar calculations are provided for 
39 K and 133 Cs, as seen in Fig. [8] and [9] respectively. We 
have shown that the MDDI can have substantial impact 
on the shape and stability for each of these alkali BECs. 



3. Other Alkalis 

While some of the other alkalis have potential for ob- 
serving the effects of MDDI, the effects arc usually small. 
For example, in 23 Na and 85 Rb BECs with atom number 
and trapping frequencies similar to current experiments 
(see Table [TTj) , a calculation including MDDI makes a five 
to ten percent difference in TOF aspect ratio from a cal- 
culation where MDDI are not included. Thus, though 
small, the MDDI effect lies within the bounds of possible 
experimental observations. For 41 K and 87 Rb, however, 
the Feshbach resonance is too narrow to provide the pre- 
cision control of a s to carry out the type of experiments 
proposed here, and the MDDI effect is less than a one 
percent perturbation on the aspect ratio. 



FIG. 4. The aspect ratio, (a), the axial length, (b), and the 
radial length, (c), for a 52 Cr BEC evolving in trap with an 
initial radial length twice the static value and the initial axial 
length half the static value, (d) and (e) show the free expan- 
sion for the BEC if it were released at times from the trap 
when the aspect ratio is at a minimum and maximum values 
respectively. The release points are indicated with arrows in 
(a). The parameters used in (a-e) are f z = f r = 700Hz, 
N = 2 x 10 4 , and a s = 14a . 



In Fig. [7j (a) shows the effect of A and a s on the sta- 
bility of a BEC with MDDI. The line for each simulated 
BEC atom numbers N indicates the boundary between 
the stable (above) and unstable (below) regimes. Similar 
to the 52 Cr case (Fig. [2}, the effect of MDDI is to stabi- 
lize a BEC in a pancake trap, and to destabilize a BEC 
in a cigar trap, (b) shows k over a range of A, with the 
collapsed regime indicated by the dotted faint line, (c) 
shows plots of k after release from a trap. The in-trap 
starts with a s = O.Olao (in the stable regime), and upon 
release (at t — in the simulation) a s is tuned to O.OOlao 
to induce collapse in free expansion. Such rapid modula- 
tion of a s has already been performed experimentally in 



IV. CONCLUSION 

We showed that the variational method provides a use- 
ful and simple tool to simulate the effects of MDDI in 
BECs and presented various results for 52 Cr and the al- 
kalis. For example, examining the the aspect ratio of 
freely expanding BECs should be sufficient for detecting 
the effects of MDDI in many species, and we suggest the 
investigation of 7 Li, 39 K , and 133 Cs as favorable to de- 
tect such effects. We mention that future investigation of 
MDDI among non-alkali species looks promising as well. 
The achievement of the magnetic-optical cooling of 168 Er 
(/x = 7fj, B ) 41j and BEC of 164 Dy (ji = 10/i B ) 9J is quite 
exciting, as the add for 168 Er and 164 Dy are 66.3 ao and 
131. 5ao, much greater than even that of 52 Cr. The abil- 
ity to tune a s by Feshbach resonance could make these 
species unparalleled for the observation of strong MDDI 
effects. The method presented here is applicable to both 
species. To close, we highlight that the effects of MDDI 
on the BEC shape provide a clear and intuitive picture 
of MDDI in BECs. The examination of the BEC aspect 
ratio, for example, can be used as a sensitive measure- 
ment of the a s value in situ, and may prove a helpful 
calibration method for future studies of other — perhaps 
more exotic — MDDI effects. We also note that while our 
discussion is limited to magnetic dipole-dipole interac- 
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FIG. 5. Simulation of a J Cr BEC frequency spectra of the in-trap oscillations. Frequency spectra are obtained by taking the 
Fourier transform of the in-trap oscillations of the aspect ratio k (a,d), radial length q r (b,e), and the axial length q z (c,f) 
that occur after changing the axial size to q a ,i from its in-trap equilibrium value, q a ,o- The colorbar gives the amplitude of the 
FFT for each plot. Note the shifts in both amplitude and frequency of oscillations between the simulations including (a-c) and 
not including (d-f) MDDI. The parameters assumed a spherical trap with f r = f a = 700Hz, N=2 x 10 4 , and a s = 14ao. 



tions in this paper, the variational method we present is 
general for all dipolar BECs and may be employed in cal- 
culating effects of electric dipole-dipole interactions (e.g. 
polar molecular BECs) as well [18] . 



Appendix A: Supplemental Experimental 
Parameters 

Typical atom numbers and trap frequencies currently 
employed in experimental studies of alkali BECs from 
the literature are listed below to show the experimental 
feasibility of the simulations provided here (Table [TI| . 



Species 


fr (Hz) 


/- (Hz) 


Num c 


f Atoms 


Rcf 


7 Li 


193 


3 


3 x 


10 s 


1 


23 Na 


1500 


150 


3 x 


10 s 


42 


39 K 


65 to 74 


92 


3 x 


10 4 


43 


41 K 


325 


15 


3 x 


10 5 


44 


85 Rb 


17 


6.8 


1 X 


10 s 


45 


87 Rb 


930 


11 


3.6 


x 10 6 


30 


133 Cs 


14 


14 


1.6 


x 10 4 


46 


52 Cr 


600 


370 


3 x 


10 4 





TABLE II. An example list of representative parameters for 
alkali BECs from the literature. 



Appendix B: Calculating the Magnetic Moment 

A key parameter for the variational simulation is the 
value of the magnetic dipole moment, as a^d °c fJ. 2 . In low 
fields, the magnetic moment is found from the Zeeman 
effect, where 

AE\ Fmp> = i-i B 9Fm F B z (Bl) 

In such a case, /i = \iB9F m F- For higher fields, however, 
the Breit-Rabi formula is used for the ground states of 
alkali atoms [47] : 

-AE hfs AE hfs ( 2m F 2 \ 1 

E = — —+q I a B m F B± — 1 H — x + x 

2(21 + 1) yP 2 ^ 1 + 1/2 J 

(B2) 

where m F = mj ± 1/2, AEhfs = Ahf s (I + 1/2) is the 
hyperfme splitting, x — , and g = gj — gi. 

The magnetic dipole moment is simply the derivative 
of the energy with respect to the magnetic field, resulting 
in 

2wtF i x 

V(B) = giVBmp ± — ^±1 ~yj^9^b (B3) 




FIG. 6. The in-trap aspect ratio for a BEC near a s = 0. The 
colored, vertical bars indicate the minimum scattering length 
currently achievable in typical experiments, a 3im i n . Calcu- 
lation for each species assumes a cigar trap with A = 100 
where f z = 20 Hz, f r = 2000 Hz, and 10 6 atoms in the 
BEC. A thick dot at the end of each solid curve indicates the 
^thrcshoid^ wn j cn depends both on the mass and the magnetic 
moment of the atom. The dotted line, shown only for 133 Cs 
as a demonstrative example, indicates the numerically unsta- 
ble solutions for a B < a * rcshold anc i j s US ed to determine the 
threshold a s . Collapse threshold values from the smallest to 
largest are (0.01, 0.02, 0.20, 0.24, 0.54, 13.50)xa for ( 7 Li, 
23 Na, 85 Rb, 39 K, 133 Cs, 52 Cr) respectively. 
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O-bg (O-O 


) Calculated ^ cross /^B 


Rof 


7 Li 


3/2 


401.752 


-0.001182 


|1, 


+ 1 > 


736.8 


-192.3 


-24.5 


0.84 


m 


23 Na 
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-0.000805 


jl, 


+ 1 > 


907 
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63 


0.50 


ED 


39 K 
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1431148] 
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TABLE III. Various parameters for the atomic species discussed in this paper. Values for gi are from Ref. [52] ■ Values for A^fs 
are from [47) ■ We use gj — 2.0023193043622(15) from |53| . Other values are from references listed in the table and [3]. * from 
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FIG. 7. Effect of MDDI in a 7 Li BEC. (a) the calculated 
threshold a s between stable and unstable regimes for a range 
of A and three representative atom numbers, (b) the BEC 
aspect ratio over a range of A, with the collapse regime in- 
dicated by the faint dotte d lin e, for which there is no stable 
solution found for Eqns. 
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(c) Evolution of the 



and 4b 

aspect ratio of the BEC in TOF expansion after release from 
a trap. In the trap a s — O.lao (where the BEC is stable), 
and upon release a„ is tuned to O.OOlao to induce collapse in 
free expansion, (d) the in-trap evolution of a BEC similar 
to (c) with a s tuned to O.OOlao at t = to induce in-trap 
dipolar collapse. Parameters used in the simulation: (a) and 
(b) favg = 700 Hz; (b), N = 5 X 10 6 and a s = 0.001a ; (c) 
and (d), f z = 200 Hz and f r = 2000 Hz with N = 5 x 10 6 . 



FIG. 8. Effect of MDDI in a 39 K BEC. (a) the calculated 
threshold a s between stable and unstable regimes for a range 
of A and three representative atom numbers, (b) the BEC 
aspect ratio over a range of A, with the collapse regime in- 
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solution found for Eqns. 4a and 4b (c) Evolution of the 



and 4b 

aspect ratio of the BEC in TOF expansion after release from 
a trap. In the trap a„ = 0.5ao (where the BEC is stable), 
and upon release a 3 is tuned to 0.05ao to induce collapse in 
free expansion, (d) the in-trap evolution of a BEC similar to 

(c) with a s tuned to 0.05ao at t = to induce in-trap dipolar 
collapse. Parameters used in the simulation: (a) and (b), 
favg = 700 Hz; (b), N = 5 X 10 5 and a s = 0.05a ; (c) and 

(d) , f z = 200 Hz and f r = 2000 Hz with N = 5 x 10 5 . 
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